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ABSTRACT 

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal 
enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primor- 
dial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption 
lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive 
the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation 
of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial 
temperature > 10 4 K. We explore the unusual circumburst environment inside the systems that hosted 
Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For 
minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) 
with values of about a few cm" 3 . In more massive halos, corresponding to the first galaxies, the 
density may be larger, n > 100 cm -3 . The resulting afterglow fluxes are weakly dependent on redshift 
at a fixed observed time, and may be detectable with the James Webb Space Telescope (JWST) and 
Very Large Array (VLA) in the near-IR and radio wavebands, respectively, out to redshift z > 20. We 
predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak 
fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected 
to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the 
early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first 
galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages 
(C II, O I, Si II and Fe II). We show that if the afterglow can be observed sufficiently early, analysis of 
the metal lines may distinguish whether the first heavy elements were produced in a pair-instability 
supernova (PISN), or a core-collapse (Type II) SN, thus constraining the initial mass function of the 
first stars. 

Subject headings: cosmology: observations - cosmology: theory - galaxies: high-redshift - gamma 
rays: bursts - quasars: absorption lines 



1. INTRODUCTION 

One of the key goals in modern cosmology is to study 
the formation of the first stars and galaxies at the end 
of the cosmic dark ages, and how they shaped the subse- 
quent evolution of the universe (Barkana & Loeb 2001; 
Bromm & Larson 2004; Ciardi & Ferrara 2005; Bromm et 
al. 2009; Loeb 2010). The first, so-called Population III 
(Pop III) stars are predicted to have formed at z > 20 
in minihalos with virial mass ~ 10 6 M© and tem- 
peratures T v ir < 10 4 K (Haiman et al. 1996; Tegmark et 
al. 1997; Yoshida et al. 2003). The Pop III initial mass 
function (IMF) is thought to be top-heavy (Bromm et 
al. 1999, 2002; Abel et al. 2002), possibly extending to 
> 100 M@, but recent simulations indicate that the 
Pop III IMF may be quite broad, also including a fraction 
of lower-mass stars (Stacy et al. 2010; Clark et al. 2011b; 
Greif et al. 2011, 2012). The first bona-fide galaxies are 
expected to have formed at a later stage in hierarchi- 
cal structure formation (Bromm & Yoshida 2011), when 
~ 10 8 M Q halos assembled at z > 10 via the merging 
of progenitor minihalos (Wise & Abel 2007; Greif et al. 



2008). These systems are often termed 'atomic cooling 
halos', because their virial temperature, T v ; r > 10 4 K, 
exceeds the threshold to enable efficient cooling via lines 
of atomic hydrogen (Oh & Haiman 2002). Direct obser- 
vations of the first galaxies at redshifts z > 10 have so far 
been out of reach. In the coming decade, the James Webb 
Space Telescope (JWST) promises to directly probe this 
critical period (Gardner et al. 2006). The detection 
of metal absorption lines in the afterglow spectrum of 
high-redshift gamma-ray bursts (GRBs), imprinted by 
enriched gas in the first galaxies, offers an unusual op- 
portunity to study the physical conditions inside them. 
We may thus be able to derive constraints on the temper- 
ature, metallicity, ionization state, and kinematics in the 
interstellar medium (ISM) of high-redshift galaxies, and 
in the surrounding intergalactic medium (IGM). It is en- 
couraging that such diagnostics can already be obtained 
for bursts at somewhat lower redshifts. An example is 
GRB 081008, where high-resolution spectroscopy with 
the VLT has probed the ISM of a host galaxy at z ~ 2 
(D'Elia et al. 2011). 
Long-duration GRBs have been shown to be associ- 
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ated with the death of massive stars (Stanek et al. 2003; 
Hjorth et al. 2003; Woosley & Bloom 2006). Their high 
luminosities make them detectable out to the edge of the 
visible universe (Ciardi & Loeb 2000; Lamb & Reichart 
2000; Bromm & Loeb 2002, 2006; Gou et al. 2004; Inoue 
et al. 2007; Meszaros & Rees 2010), with the current 
record held by GRB 090429B at z ~ 9.4 (Cucchiara et 
al. 2011). GRBs provide ideal probes of the high- redshift 
universe, including the star formation rate (Totani 1997; 
Wijers et al. 1998; Porciani & Madau 2001; Chary et 
al. 2007; Yiiksel et al. 2008; Wang & Dai 2009, 2011; 
Elliott et al. 2012), reionization (Gallerani et al. 2008), 
dark energy (Dai et al. 2004; Wang et al. 2011), and the 
IGM metal enrichment (Barkana & Loeb 2004; Totani et 
al. 2006; Toma et al. 2011; Bromm & Loeb 2012). The 
leading contender for the central engine of long-duration 
GRBs is the collapsar model (Woosley 1993; MacFadyen 
et al. 2001). Because of their predicted high characteris- 
tic mass, a significant fraction of Pop III stars might end 
their lives as a black hole, potentially leading to a large 
number of high-redshift GRBs. Thus, Pop III stars are 
viable progenitors of long-duration GRBs, triggered by 
the collapsar mechanism, as long as they can lose their 
outer envelope and retain sufficient angular momentum 
in their center (Bromm & Loeb 2006; Belczynski et al. 
2007; Komissarov & Barkov 2010; Stacy et al. 2011). It 
might even be possible for Pop III collapsars to occur if 
the extended outer envelope were not lost (Suwa & Ioka 
2011). 

The history of pre-galactic metal enrichment has sev- 
eral important consequences for structure formation 
(Madau et al. 2001; Karlsson et al. 2012). An early 
phase of metal injection may qualitatively change the 
character of star formation, from a predominantly high- 
mass (Pop III) mode to a normal, low-mass dominated 
(Pop I /II) one, once the enrichment has exceeded a 
'critical metallicity' of -Z C rit ~ 1O~ 4 Z (Bromm et al. 
2001a; Schneider et al. 2002, 2006; Bromm & Loeb 2003; 
Mackey et al. 2003). The transition between these two 
modes has crucial implications, e.g., for the expected red- 
shift distribution of GRBs (Bromm & Loeb 2006; Camp- 
isi et al. 2011; de Souza et al. 2011), for reionization 
(Cen 2003; Wyithe & Loeb 2003; Furlanetto & Loeb 
2005), and for the chemical abundance patterns of low- 
metallicity stars (Qian & Wasserburg 2001; Frebel et al. 
2007, 2009; Tumlinson 2010). It is therefore important 
to explore the topology of early metal enrichment, and 
to determine when particular regions in the universe be- 
come supercritical (Tornatore et al. 2007; Maio et al. 
2010). 

Absorption lines imprinted on the spectra of bright 
background sources, such as GRBs or quasars, are one of 
the main sources of information about the physical and 
chemical properties of high-redshift systems (Oh 2002; 
Furlanetto & Loeb 2003; Oppenheimer et al. 2009). 
These lines are due mainly to absorption by neutral hy- 
drogen present in the low column-density Lya forest, and 
by metals in low-ionization stages (e.g., C II, Si II, Mg II, 
Fe II, O I) which arise in the higher column-density gas 
associated with Damped Lya Absorbers (DLAs). The 
analysis of the spectrum of distant GRB 050904 (Totani 
et al. 2006) has resulted in a wealth of detailed insight 
into the physical conditions within the host galaxy at 



z ~ 6.3, and Salvaterra et al. (2009) claimed that they 
identified two absorption lines (Si IV and Fe II), although 
at poor signal-to- noise, in the spectrum of GRB 090423, 
the most distant spectroscopically confirmed burst at 
z = 8.2 (Salvaterra et al. 2009; Tanvir et al. 2009). 
GRBs as background sources offer a number of advan- 
tages compared to traditional lighthouses such as quasars 
(Bromm & Loeb 2012). Their number density drops 
much less precipitously than quasars at z > 6 (Fan et 
al. 2006), and the absence of a strong proximity effect, 
together with the near power-law character of their spec- 
tra, renders them ideal probes of the early IGM. 

In this paper, we discuss the observational signatures 
of Pop III GRBs and study pre-galactic metal enrich- 
ment utilizing absorption lines in the spectra of high-z 
GRBs which were imprinted by the first galaxies. Re- 
cently, it has become feasible to study the formation of 
the first galaxies, including the metal enrichment from 
Pop III supernovae (SNe), with highly-resolved cosmo- 
logical simulations (Wise & Abel 2008; Greif et al. 2010). 
We here place a bright GRB into the simulation box of 
Greif et al. (2010), and derive the spectral signature as 
the afterglow light escapes from the first galaxy, thereby 
probing the partially enriched IGM in its vicinity. The 
fluxes in the near-IR and radio bands may be detectable 
by the JWST and the Very Large Array (VLA) out to 
z > 20. The structure of this paper is as follows. In 
Section 2, we derive the circumburst density profile of 
Pop III GRBs, followed by a brief description of the un- 
derlying first galaxy simulation (Section 3). We discuss 
the properties of the GRB afterglow in Section 4, and the 
metal absorption line diagnostics in Section 5, followed 
by our conclusions. 

2. GRB DENSITY PROFILE 

Two different environments are currently discussed as 
Pop III star formation sites, minihalos and atomic cool- 
ing halos (Bromm et al. 2009). The latter mode is often 
termed Pop III. 2 to indicate physical conditions, such 
as a higher degree of ionization, that differs from the 
canonical minihalo case, possibly resulting in somewhat 
lower masses (McKee & Tan 2008). For normal GRBs, 
either a constant number-density profile, or a power-law 
dependence, as might be expected in the stellar wind 
from the progenitor, is discussed in the literature (Bloom 
2011). The GRB afterglow emission sensitively depends 
on the circumburst density (Ciardi & Loeb 2000; Gou et 
al. 2004; Inoue et al. 2007), and this is one of the crucial 
uncertainties in making predictions for Pop III bursts. 
We next derive the likely circumburst conditions in the 
two possible Pop III host systems. 

2.1. Minihalo case 

For simplicity, we assume that Pop III GRBs are trig- 
gered in minihalos close to the cooling threshold for 
collapse, which is only weekly dependent on redshift 
(Yoshida et al. 2003). For definiteness, we choose 
M v ; r ~ 10 6 M Q . To determine the pre- burst density, we 
need to consider the build-up of the H II region around 
the central Pop III star, and how the density structure 
is modified through the strong photoionization-heating 
inside of it (Whalen et al. 2004) . The production of ion- 
izing photons strongly depends on the stellar mass, which 
in turn is determined by how the accretion flow onto the 
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growing protostar proceeds under the influence of this 
radiation field (e.g., McKee & Tan 2008; Hosokawa et 
al. 2011; Smith et al. 2011; Stacy et al. 2012). Thus, 
the assembly of Pop III stars and the development of an 
H II region around them proceed simultaneously and af- 
fect each other. The shallow potential wells in the host 
minihalos, with corresponding circular velocities of a few 
kms -1 , are unable to retain photo- ionized gas, so that 
the gas is effectively blown out of the minihalo. The 
resulting photo-evaporation has been studied with 3D 
radiative transfer calculations (Alvarez et al. 2006; Abel 
et al. 2007; Greif et al. 2009), where one massive Pop III 
star at the center of the minihalo acts as an embedded 
point source, that also take into account the hydrody- 
namic response of the photo-heated gas. 

It is possible to understand the key physics of the pho- 
toevaporation from minihalos with the self-similar solu- 
tion for a champagne flow (Shu et al. 2002). Assuming 
a p oc r~ 2 density profile, which describes the typical sit- 
uation in minihalos outside a nearly flat inner core, one 
can reformulate the spherically symmetric continuity and 
Euler equations for isothermal gas as follows: 
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where x — r/c s t is the similarity variable, and p(r,t) = 
a(x) I AirGt 2 = m,B_n(r)/X and u(r,t) — c s v{x) introduce 
the reduced density and velocity, respectively. c s is the 
sound speed of the ionized gas and X = 0.75 the hydro- 
gen mass fraction. The density globally decreases with 
time, after the ionizing source has turned on. To derive 
the immediate pre-explosion value, we will therefore set 
t = f* ~ 3 x 10 6 yr, the typical lifetime of a massive 
Pop III star (e.g., Bromm et al. 2001b). 

Different solutions are obtained depending on the ratio 
e = (c s j/c s ) 2 , where c Si i and c s are the initial and ionized 
isothermal sound speeds. We take T = 3 x 10 4 K and 
Ti = T v ir, where the virial temperature is defined as 

Tyi ^GM^ / M Y 3 (l±z \ 

2R vir k B \10 8 MqJ \ 10 J 

The intuition behind the latter identification is that 
c 2 i ~ fcelvir/mH — «virj the specific gravitational po- 
tential energy of the gas inside a halo, is a measure of 
how strongly the fluid is bottled up, before a champagne 
flow can occur. Using table 1 in Shu et al. (2002), we can 
obtain the density profiles at different redshifts as follow. 
First, the value of e can be calculated at the given red- 
shift. We then derive the corresponding reduced central 
density, olq, via interpolation. For the given ao we can 
numerically solve equations (p} and © subject to the 
inner boundary conditions: 

a = ao and v = 0, at x — 0. (4) 

To avoid the singularity near x = 0, Shu et al. (2002) 
provide the series expansion 

a = a + ^ ( \ - a ) x 2 + . . . , (5) 
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Fig. 1. — Minihalo circumburst environment. Shown is the hy- 
drogen number density vs. distance from the central Pop III star 
at the moment of its death, 3 Myr after it was born. The density 
profiles are given by the Shu solution for different redshifts, as la- 
belled. Typical circumburst densities are ~ 1cm -3 , with smaller 
values at lower redshifts. The redshift trend reflects the shallower 
minihalo potential wells at later times. 

v = l x+ h{h aa ) x3 + -- (6) 

Once the Pop III star has turned on at t — 0, the photo- 
heating commences, putting in place a pressurized inner 
bubble which leads to an outward-flowing shock at x s ^. 
The shock obeys the usual isothermal jump conditions, 
written in terms of the reduced variables as 

(v u - x sh )(v d - x sh) = 1, — = (v u - x sh ) 2 , (7) 

where the subscript u(d) indicates the value of the re- 
duced velocity and density upstream (downstream) of 
the shock. 

In Figure Q3 we show the resulting density profiles at 
the end of the Pop III progenitor's life, and therefore 
at the completion of the photoionization-heating as well. 
As can be seen, the circumstellar densities have dropped 
significantly from the high values, > 10 4 cm~ 3 , present 
prior to Pop III star formation, and are nearly uniform 
at small radii. Such a flat density profile is markedly 
different from that created by stellar winds. The ab- 
sence of any wind around a Pop III star is indeed pre- 
dicted, as a consequence of very low metallicity (Baraffe 
et al. 2001; Kudritzki 2002), and could serve as an indica- 
tor of a low-Z GRB progenitor. The post-photoheating 
density evolves with redshift, approximately according 
to n p i oc (1 + z), normalized such that n ~ lcm~ 3 at 
z ~ 20. This law differs from earlier estimates of the 
Pop III circumburst evolution, where no z-dependence 
was assumed, or a scaling with the density of the back- 
ground universe, n p ; oc (1 + z) 3 (e.g., Gou et al. 2004). 

This relation can be understood as follows. Prior to 
the onset of protostellar collapse, the baryonic density 
profile in the minihalo can be described by the Lane- 
Emdcn equation for an isothermal sphere 
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Fig. 2. — Atomic-cooling-halo circumburst environment. The 
density profiles are again calculated from the Shu solution, evalu- 
ated at t = 3 X 10 e yr after source turn-on. The collapse redshift 
is here indirectly given by the virial temparaturc, as labelled. We 
show the case of photoheating from only a single Pop III star (solid 
lines), and that from a stellar cluster (dotted lines). The result- 
ing densities again exhibit a nearly uniform inner core, but overall 
values are much higher than in the minihalo case (see Fig. 1). 

where po is the density at which the gravitational poten- 
tial <j) is set to zero. A simple analytical solution to this 
equation is the well-known singular isothermal sphere 
(SIS): psis = a 2 /(2nGr 2 ), where a 2 ~ k B T vir /m H . This 
law approximately describes the profile found in realis- 
tic simulations of first star formation (e.g., Yoshida et 
al. 2003), outside an inner core of radius rj ~ lpc, the 
Jeans length of the primordial gas (Bromm et al. 2002). 
The baryonic mass inside the shock radius, r s h = c s f* ~ 
100 pc, is about 

Mi = / 4irr 2 p SIS {z,r)dr . (9) 

Jr. j 

At the end of photoionization-heating, the mass in the 
flat inner core is about 

47r Air 
M f = Y r shP P i = y r 3 sh n pi mnX. (10) 

If we neglect any outflows at large radii, we have approx- 
imate mass conservation, such that M, ~ Mf, resulting 
in n pl cx T v - lr cx (1 + z), as before. 

Our post-photoheating densities are consistent with 
simulation results (Kitayama et al. 2004; Alvarez et 
al. 2006; Abel et al. 2007; Greif et al. 2009, 2010). 
We also agree with Whalen et al. (2004), if we take 
into account that their halo had a somewhat lower mass, 
My- 1T ~ 5 x 10 5 M Q . Empirical support is provided 
in Chandra et al. (2010) who found that the circum- 
burst density of GRB 090423 at z = 8.2 was about 
n r~j 0.9 cm' 3 by fitting the radio, X-ray and infrared 
afterglow, although this burst almost certainly did not 
originate in a minihalo environment. 

2.2. Atomic cooling halo case 

Star formation and radiative feedback inside atomic 
cooling halos is significantly less well-understood than 



the minihalo case (e.g., Johnson et al. 2009; Safranek- 
Shrader et al. 2010, 2012), and there are no comparable 
high-resolution simulations yet. What is the character of 
star formation, in terms of IMF and stellar multiplicity 
(Clark et al. 2011a)? Given current uncertainty, we here 
again use the formalism of the Shu solution as developed 
above. We further assume for simplicity that either a 
single Pop III star forms, or a small stellar cluster, with 
properties that would again heat the surrounding H II 
region to T ~ 3 x 10 4 K. We stress that the numbers 
derived here only provide us with rough guesses at best, 
and the hope is that simulations will eventually become 
available to firm them up. 

If a cluster of Pop III stars forms, the time where pho- 
toheating stops will be prolonged. We thus need to eval- 
uate p(r,t) = a(x) I AirGt 2 at t = truster — £*+^sf, where 
tsF is the timescale over which star formation is ongoing, 
and t* again is the lifetime of a single Pop III star. We 
can estimate tsF ~ t*, because disruptive feedback ef- 
fects will tend to terminate star formation when the first 
member stars die. The resulting densities at the end of 
photoheating are shown in Figure [2j Similar to the mini- 
halo case, densities are nearly constant at small radii, but 
overall values are much higher. The reason is that here 
the potential wells are deeper, so that photoheated gas 
can more easily be retained. Typical circumburst den- 
sities are n p i ~ 100 cm -3 . Because such large densities 
would correspond to high afterglow fluxes, Pop III GRBs 
originating in atomic cooling halos may be extremely 
bright, rendering them visible out to very high redshifts. 
For some of the most distant bursts known, such large cir- 
cumburst densities have indeed been inferred. E.g., Gou 
et al. (2007) have argued that fitting the GRB 050904 
afterglow at z = 6.3 requires a circumburst density of 
a few 100 cm~ 3 . Similar values have been suggested for 
GRB 080913 at z = 6.7 (Greiner et al. 2009; Zhang et 
al. 2009). 

3. MODELING PRE-GALACTIC ENRICHMENT 

In order to predict afterglow spectra from Pop III 
GRBs, we consider explosions that are embedded in 
a realistic cosmological setting. Specifically, we use 
the first galaxy simulation carried out by Greif et al. 
(2010), where the assembly of an atomic cooling halo 
was tracked, resolving all prior Pop III star formation 
in the progenitor minihalos, together with the radiative 
feedback from H II regions around those stars. In addi- 
tion, the simulation allowed one of the Pop III progenitor 
stars to explode as an energetic SN, explicitly following 
the transport and mixing of the ejected metals into the 
IGM. We here summarize a few key aspects of the sim- 
ulation, and refer the reader to Greif et al. (2010) for 
further details. 

The simulation is carried out in a cosmological box of 
size 1 Mpc (comoving), and is initialized at z = 99 ac- 
cording to a concordance A cold dark matter (ACDM) 
model with matter density Q m = 1 — Q\ = 0.3, 
baryon density Sib = 0.04, Hubble parameter h = 
H / (100 km s _1 Mpc -1 ) = 0.7, spectral index n s = 1.0, 
and normalization as = 0.9 (Spergel et al. 2003). Den- 
sity and velocity perturbations are imprinted at recombi- 
nation with a Gaussian distribution, and they are propa- 
gated to z = 99, when the simulation is started, by apply- 
ing the Zeldovich approximation. Within the highest res- 
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Fig. 3. — Possible explosion sites for high-z GRBs. Shown are the hydrogen number density and metallicity contours during the assembly 
of a first galaxy, averaged along the line of sight within the central ~ 100 kpc (comoving), at two different rcdshifts. Panel (a): z ~ 25, 
briefly after the first Pop III SN exploded. At this time, most of the metals reside in the 1GM. Panel (b): z ~ 16.4, closer to the virialization 
of the atomic cooling halo. Now, metals are being re-assembled into the growing potential well of the first galaxy. The topology of metal 
enrichment is highly inhomogeneous, with pockets of highly enriched material embedded in regions with a largely primordial composition. 
These plots are derived from the simulation described in Greif et al. (2010). 



olution region of our nested, zoomed-in initial conditions, 
the DM and gas particle masses are 33M Q and 5M , re- 
spectively, The baryonic mass resolution, roughly given 
by the mass contained inside an SPH kernel, is ~ 400M©, 
close to the relevant Jeans mass in the primordial star 
forming gas. To capture the chemical evolution of the 
gas, the simulation follows the abundances of H, H + , H~, 
H 2 , Hj, He, He + , He ++ , and e~, as well as the five deu- 
terium species D, D + , D — , HD and HD + . All relevant 
cooling mechanisms are included, i.e., H and He atomic 
line cooling, bremsstrahlung, inverse Compton scatter- 
ing, and collisional excitation cooling via H 2 and HD 
(Johnson & Bromm 2006). 



In Figure [3l we show the hydrogen number density and 
metallicity averaged along the line of sight within the 
central ~ 100 kpc (comoving) at two different output 
times, from z ~ 25, briefly after the first star-forming 
minihalo has collapsed, to z — 16.4, closer to the viri- 
alization of the first galaxy. The distribution of metals 
produced by the first SN explosion is highly inhomoge- 
neous, and the metallicity can reach up to Z ~ 10~ 2 ' 5 Zq, 
which is already larger than any of the values discussed 
for the critical metallicity, i? C rit < lO^Z© (e.g., Bromm 
et al. 2001a; Bromm & Loeb 2003; Frebel et al. 2007), in- 
cluding those predicted for a dust-driven Pop Ill/Pop II 
transition, where 10~ 6 Z & < Z CI a < 1O~ 5 Z (Omukai et 
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Fig. 4. — Environment for high-2 GRBs. Shown are tempera- 
ture contours along the line of sight within the central ~ 100 kpc 
(comoving) at z = 16.4, corresponding to panel (b) in Fig. [3] The 
center of the growing galaxy is highly turbulent, with pockets of 
cold gas embedded in a hotter, diffuse medium. 

al. 2005; Schneider et al. 2006). The implication is that 
both Pop III and Pop I /II stars will form during the as- 
sembly of the first galaxies (Johnson et al. 2008; Maio 
et al. 2010), in turn giving rise to the simultaneous oc- 
currence of Pop III and normal GRBs at a given redshift 
(Bromm & Loeb 2006; de Souza et al. 2011). Here, we 
will focus on the case of a Pop III burst exploding in one 
of the (still metal-free) first galaxy progenitor minihalos 
at z ~ 16.4 (see Section 5). In Figure |H we show the 
temperature distribution in the vicinity of the explosion 
site, allowing us to calculate (thermal) line broadening 
in a realistic fashion. 

4. POP III AFTERGLOW EMISSION 

4.1. Basic Physics 

Our modeling of the GRB broad-band afterglow emis- 
sion follows standard prescriptions, and we here only 
briefly summarize the basic physics involved (for general 
reviews, see Meszaros 2006; Zhang 2007). We specifi- 
cally consider a relativistic shell ejected from the pro- 
genitor of a GRB (e.g., Kobayashi 2000). The shell has 
a rest mass of Mq, an observed thickness of Ao, and 
exhibits an initial Lorentz factor of Tq = Ei SO /Moc 2 , 
where i?i so is the isotropic-equivalent explosion energy. 
The interaction between the shell and the surrounding 
medium with a density of no produces a forward and 
a reverse shock. The shocks, in turn, accelerate elec- 
trons to high energies, distributed according to a power- 
law: dN{ r ) e )/d r y e cx j~ p for 7 e > 7 TO , where j e is the 
internal Lorentz factor for the electrons. The minimum 
Lorentz factor is set by assuming that a fraction e e of the 
post-shock energy is transferred to the electrons, result- 
ing in: 7 m = e e j(m p /m e )(p - 2)/{p - 1), where 7 ~ T 
is the Lorentz factor of the shocked fluid (e.g., Blandford 
& McKee 1976; Sari et al. 1998). If one also assumes 
that a fraction eg of the post-shock energy is deposited 




Log t (s) 



Fig. 5.— Flux density at v = 6.3 X 10 13 Hz (M band) as a func- 
tion of observed time. We show the emission from the forward 
shock (dashed line), the reverse shock (solid red line), and their 
combination (solid black line). We use the parameters given in the 
text, for a GRB exploding in a minihalo with circumburst den- 
sity no = 0.5 cm -3 at z = 10. It can be seen that the reverse 
shock dominates during the first few hours, reaching a brightness 
in excess of ~ 1 mjy. 

into random magnetic fields, the gyration of the shock- 
accelerated electrons around the magnetic field will gen- 
erate synchrotron radiation, thus giving rise to the af- 
terglow emission. There is then a second characteris- 
tic Lorentz factor, j c , determining the threshold below 
which electrons do not lose a significant fraction of their 
energy to radiation. Equating the synchrotron cooling 
time with the expansion time of the shell, one finds: 
7c m e c 2 = P(7cKxp, where P( 7e ) = (4/3)a T c7 e 2 (B 2 /87r) 
is the synchrotron radiation power (Rybicki & Lightman 
1979), and t exp the expansion time, measured in the co- 
moving fluid frame (Sari et al. 1998; Panaitescu & Ku- 
mar 2000). Therefore, the cooling Lorentz factor can be 
written as: j c = 6irm e c/(aTB 2 t cxp ). 

For the resulting spectrum, there are now two cases: In 
the fast cooling case ("f rn > j c ), the flux can be calculated 
from 



F u = F„ 



(^) 1/3 (^) 2 

(Ji.U/3 



V < v a 
V a < V < v c 
V c < V < v m 
V m < V, 



whereas in the opposite, slow-cooling (^f m < 7 C ), case, 
one has 



F„ = F u 



(f^) 1/3 (^) 2 

(_^)l/3 

(J^)-(p-1)/2 

(|£.)-(p-1)/2(^)-J>/2 



V < V a 

V a < V < V „ 
V m < V < v c 

V r < V. 



The break frequencies correspond to the electron 
Lorentz factors introduced above. The frequency v a de- 
scribes the onset of synchrotron self-absorption, which 
is important at low (radio) frequencies. The afterglow 
emission exhibits a different character before and af- 
ter the reverse shock (RS) has crossed the shell. This 
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happens at a time f® ~ (1 + z)A /2c = 16.7s (1 + 
z)A /10 12 cm, as measured in the observer frame. Be- 
fore RS crossing (t < i®), we have for the forward shock 
07"): 

(1 i \-l 2 1/2 a -3/2^1/2,-1 

/-, , \-l -3/2 -1.-1/2^-1/2,-1 

i/ Cl/ oc (1 + z) e B n A £ iao ' * , 



y a,/« (1 + z)~ 



-1^6/5^11/10^-1/2 p 7/10 + -2 



& iso 1 



-F^max,/ oc (1 + z)e 1 J 2 n 1 /2 E iso d L 2 t, 



(11) 

where dj, is the luminosity distance. Here and in the 
following we equate the pre-shock density with that es- 
tablished by photoionization-heating, at the time of the 
Pop III star's death: no ~ n p \ (see Section 2). For the 
reverse shock ("r"), we have 

v m ,rcc (1 + zy 1 e 2 e e 1 J 2 Tln 1 /2 , 
i/^cx (l + zrh B 3/ WA- 1/2 Er s ^ 2 t-\ 
u a , r oc (1 + zrhi^^T-^A^E^h-"/ 5 , 

F v , max , r <x (l + z^nl^A-^E^dlV/ 2 . (12) 

After RS crossing, we similarly have 



VcJ = Vcj{t®) — 

. Iffi 



-3/2 



-1/2 



max,/ 



' uiax,/(*ffi) I — 



(13) 



for the forward shock. For the reverse shock, we have 



t 



v c,r — v c,r{t®) ( — 
tat 



V a ,r — Va,r(t@) I — 
tm 



-3/2 



-3/2 



-1/2 



tf 



(14) 



The exact expressions, and further details, can be found 
in Kobayashi (2000) and Wu et al. (2003). Note that the 
latter study assumes the presence of a circumburst wind 
profile, whereas we here consider the constant-density 
(ISM) case, similar to Kobayashi (2000). The basic pa- 
rameters of the afterglow emission region are: Tq — 300, 



10" erg, A 



10 



12 



0.3, e B = 0.1, 



and p — 2.5. We have verified that our afterglow pre- 
dictions are consistent with those in Gou et al. (2004), 
and Inoue et al. (2007), as these papers focus on Pop III 
bursts at high redshifts as well. As an example, in Fig- 
ure [5j we show the M-band light curve, corresponding to 
v = 6.3 x 10 13 Hz. It can be seen that the reverse shock 



t=6min 



t=1hour 



o -o -o 



o -o -o -o-o-o- 



t=1day 



o o-o-o ooooo 



10 12 
redshift z 



Fig. 6.— Observed flux at v = 1.36 X 10 14 Hz (K band) as a 
function of redshift. The lines indicate the flux at different observed 
times, as labelled. Here, we have assumed that the Pop III GRB 
is triggered in a minihalo environment, as modelled in Section 2.1. 
We mark the K-band sensitivity of the NIRSpec instrument on 
board the JWST with a horizontal line. The sharp cut-off at z ~ 17 
is due to hya absorption in the IGM. 
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Fig. 7. — Observed radio flux (in mjy) at v = 5 GHz as a function 
of redshift, again for a minihalo burst environment. We employ the 
same convention for the lines and symbols as in Fig. [6] but we now 
mark the corresponding sensitivity {horizontal line) for the EVLA. 
As is evident, the radio flux increases with observed time, and is 
well within reach of detection. 

emission dominates during the first few hours, whereas 
the forward shock component takes over at later times. 

4.2. Redshift Dependence 

We here briefly summarize the key aspects of after- 
glows triggered by Pop III stars. Figure [B] shows the 
observed flux at v = 1.36 x 10 14 Hz as a function of red- 
shift for Pop III GRBs in a minihalo. According to Sec- 
tion 2.1, circumburst densities can then be written as 
n ~ lcm~ 3 (l + z)/20. The lines with filled dots, black 
triangles and open dots correspond to an observed time 
of 6 minutes, 1 hour and 1 day respectively. The straight 
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line marks the K-band sensitivity for the near-infrared 
spectrograph (NIRSpec) on JWST, estimated for a res- 
olution of R = 1000, a signal-to-noise ratio of S/N= 10, 
and an integration time of 1 hour (Gardner et al. 2006). 
Once the observed frequency corresponds to an emitted 
frequency above the Lya resonance, v a = 2.47 x 10 15 Hz, 
which occurs at z ~ 17 for the K-band, all flux will be 
completely absorbed by the intervening, still largely neu- 
tral IGM. As can be seen, in the K band the JWST will 
be able to detect GRBs, and to conduct spectroscopy on 
their afterglows, out to z ~ 16 even after 1 day. In the 
M band, the redshift horizon is extended further still, 
to z ~ 35. At these frequencies, the afterglow flux is 
only weakly dependent on redshift. The main reason is 
that cosmic time dilation implies that a given observed 
time after the trigger corresponds to successively ear- 
lier emission times, where intrinsic afterglow luminosities 
are much brighter (Ciardi & Loeb 2000; Bromm & Loeb 
2012). Secondly, circumburst densities modestly increase 
with redshift (see Section 2). 

In Figure [3 we show the flux at v = 5 GHz, again 
after t = 6 minutes, 1 hour and 1 day. We calculate the 
sensitivity of the VLA, using the following expression 
(Ioka & Meszaros 2005) 



S/N • 2k B % 



ays 



A cS ^/2t int Ai/ 



-1/2 



lday/ 



Av 



50 MHz 



-1/2 



(15) 



where U nt is the integration time, Av the bandwith, and 
A cS /T sys - 2 x 10 6 cm 2 K -1 . The bandwidth of the Ex- 
panded VLA (EVLAjE can be up to 8 GHz. For a 2h 
integration in the 8 GHz band, the EVLA can thus reach 
a sensitivity of up to 6.4 /xJy, which is shown as a horizon- 
tal line in Figure [7] We can see that the radio afterglows 
are easily within reach for the EVLA, up to z ~ 40 (e.g., 
Chandra et al. 2010). 

4.2.1. Peak fluxes and wavelengths 

According to hierarchical cosmic structure formation, 
the characteristic halo mass varies as a function of red- 
shift. Thus, the typical observed peak flux and wave- 
length of Pop III GRBs are also a function of redshift. 
In the following, we derive ballpark estimates for these 
key quantities by considering 2a overdensities (e.g., fig- 
ure 4 in Clarke & Bromm 2003). This provides us with 
the typical halo mass, M v i r , at a given redshift. In the 
2a case, halo masses become less than 1O 5 M at redshift 
z > 15. Since this has dropped below the threshold mass 
for H 2 cooling (e.g., Yoshida et al. 2003), Pop III stars 
could not form in such low-mass halos. For higher-sigma 
peaks, however, super-critical halos still exist at earlier 
times. Therefore, for simplicity we adopt M v ; r = 10 5 M Q 
at z > 15. After the characteristic halo mass is thus 
determined, the virial temperature is fixed as well, and 
the ratio e = (c s .i/c s ) 2 is known. The Shu solution then 
gives the circumburst densities at a given redshift (see 
Section 2), and GRB afterglow spectra can finally be 
calculated. Figure [8] shows the observed peak flux at ob- 
served time t = 6 minutes, 1 hour and 1 day at z < 30. 

1 http://www.aoc.nrao.edu/evla/ 
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Fig. 8. — The observed maximum flux density as a function of 
redshift at different times, t = 1 day, 1 hour and 6 minutes from 
bottom to top. We assume that the Pop III GRB has exploded in 
a minihalo. As can be seen, afterglow fluxes can reach substantial 
levels early on, and still remain at ~ mjy a day after the trigger, 
out to very high redshifts. 
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Fig. 9. — The observed maximum wavelength as a function of 
redshift at different times, as labelled, for a GRB exploding in a 
minihalo. The emission maxima range from mm to mid-IR wave- 
lengths. 

The peak flux varies from milli-Jansky to sub-Jansky. 
Figure |9] gives the observational peak wavelength at the 
same observed times for the same redshift range. The 
observed peak wavelengths vary from millimeter to mid- 
infrared bands. Early afterglows are best observed with 
facilities such as the JWST, or possibly adaptive-optics 
supported next-generation, extremely-large telescopes on 
the ground. At late times, radio facilities provide the op- 
timal chance for detection. Among them is ALMA, which 
will allow to probe the afterglow in the mm and sub-mm 
regime, where the peak of the emission is located about 
1 day after the explosion for z ~ 10 sources. 

5. METAL ABSORPTION LINES 

As we have seen above, Pop III GRBs can in principle 
be detected out to very high redshifts. Their afterglow 
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spectra provide us with a unique probe into the state 
of the early IGM, including its degree of ionization (e.g., 
Totani et al. 2006) and metal enrichment. The latter will 
be encoded in absorption lines imprinted on the smooth 
afterglow emission with a signature that depends on the 
nature of the Pop III SNe. Here, we will consider for 
simplicity that prior to the GRB only one nearby SN ex- 
ploded beforehand, dispersing its complement of heavy 
elements into the otherwise pristine IGM. Such a scenario 
is consistent with the clustered nature of Pop III star 
formation, as is expected for the high-cr peak host sys- 
tems (e.g., Greif & Bromm 2006). Specifically, we here 
explore nucleosynthetic yields for Type II core-collapse 
SNe, and for pair- instability supernovae (PISNe). In Ta- 
ble 1, we summarize the average yields per SN, Yx, for 
select elements that are expected to dominate the spec- 
tral signature. For Type II SNe, we use the yields of 
Woosley & Weaver (1995). These depend on the explo- 
sion energy, which introduces uncertainties of ~ 25%. 
For the PISN scenario, we utilize the yields of Heger & 
Woosley (2002; 2010). Table 1 also lists the properties of 
several key transitions, including their rest-frame wave- 
length, A;, and their oscillator strength, f osc ,i- 

We employ a simple model for the ionization structure 
of the metal-enriched gas. If the hydrogen in this region 
is substantially neutral, metals will reside in states typ- 
ical of H I regions in the Milky Way (C II, O I, Si II, 
and Fe II for the elements of interest), because photons 
able to further ionize these elements will be absorbed by 
H I (Furlanetto & Loeb 2003). The ionization poten- 
tial of O I is 13.62 eV, and it remains locked in charge 
exchange equilibrium with H I (Oh 2002). A modest 
extragalactic radiation background at frequencies below 
the Lyman limit could easily maintain these metals in a 
singly- ionized state (e.g., Furlanetto & Loeb 2003). We 
thus assume that C, O, Si, and Fe are predominantly in 
the ionization states typical for Galactic H I clouds. 

5.1. Cosmological Radiative Transfer 

Employing comoving coordinates, the equation of cos- 
mological radiative transfer is (e.g., Abel et al. 1999): 



1 dl v , n • VI V H(t) , 8J V , 
c at a c av 



- h T 



(16) 



where l v is the intensity of the radiation field, h a unit 
vector along the direction of the ray, H(t) = a/a the 
Hubble constani0, and a = a(t + dt)/a(t) the ratio of 
cosmic scale factors separated by dt. Here j v and fc„ 
denote the emission and absorption coefficient, respec- 
tively. Equation (fT6| can be simplified if the effects 
of cosmic expansion are negligible. This will be the 
case if the light-crossing time, <n g ht = L/c, over the re- 
gion of interest, here the (physical) size of our simula- 
tion box, L = lMpc/(l + z), is small compared to the 
Hubble time, t H ~ 5 x 10 8 yr[(l + z)/10]" 3 / 2 . Since 
flight <C £h, the transfer equation reverts to its non- 
cosmological form: 



ldl v 



fl ■ VI V = j V — kyly 



(17) 



c dt 

2 We here employ the same cosmological parameters as in the 
underlying simulation (see Sec. 3). 



We further neglect any diffuse (re-)emission of radiation, 
such that j v = 0, and assume near steady-state condi- 
tions (d/dt ~ 0), resulting in the static transfer equation: 



(18) 



If ly(0) is the intensity of the central afterglow, approx- 
imately treated as a point source, we have the standard 
solution: l v = 1^(0) exp(— t„), where r v is the optical 
depth along the line-of-sight (LOS). 

We can obtain the optical depth as follows. For a 
(physical) differential distance element, dl, along the 
given LOS, we have: dr v = a u nxdl, where 



ov = \fne 2 /{m e c)f OSCfi H{u,x)/Au D 



(19) 



is the frequency dependent absorption cross-section, and 
nx the number density of species X (e.g., O I, C II, Fe II, 
Si II, H I). In the above equation, Afrj = bvi/c is the 
Doppler widtlfl with parameter b = ^2k^T '/mg. The 
Voigt function, H, is given by 



H(u, x) 



[x — y) 2 + u 2 



dy. 



(20) 



Here, the relative strength of natural (damping) to 
Doppler broadening is described by the parameter 



(21) 



47rAi/ D 



where Ti is the damping constant for transition i (Mor- 
ton & Smith 1973). The variables x = {v — Vi)/Av-o 
and y = v/b are the frequency difference relative to the 
line center in Doppler units, and the normalized parti- 
cle velocity, with respect to the GRB location, respec- 
tively. The total LOS value of the optical depth is then: 
T„ = J o v nxdl. We note that we can neglect any wave- 
length differences between the GRB source and the metal 
absorber, because of the very small expansion-generated 
rcdshift on the scale of our simulation box. Specifically, 
we estimate that Az = At/ (7. 5 x 10 7 yr)~ 0.005, where 
we assume z ~ 10, and At ~ tnght ~3x 10 5 yr. 

We extract the optical depths along random lines 
of sight with the ray-tracing algorithm of Greif et al. 
(2009). We choose Ng = 160 and = 320 in the spher- 
ical coordinates, which corresponds to 51200 rays. For 
each ray, we use 200 logarithmically-spaced radial bins, 
and approximate the optical depth integral with a sum 
along the LOS. The Pop III GRB will likely explode in a 
still substantially neutral IGM. Any flux shortward of the 
Lya resonance is therefore completely absorbed, and be- 
cause of the corresponding very large optical depth, there 
will be absorption in the red damping wing as well. We 
model this with the analytical formula given by Miralda- 
Escude (1998), using z rc i on = 7 and standard cosmolog- 
ical parameters (Komatsu et al. 2011). The line shape 
is not very sensitive to the choice of z rc ion (Bromm et al. 
2001b). 

5.2. Synthetic Spectra 



3 We here neglect any broadening from turbulent motions. This 
assumption is good in the minihalo progenitors, but breaks down 
in the emerging atomic cooling halo. 
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TABLE 1 

Supernova Yields and Important Transitions. 



Element 


Y x (Scalo) a 


Y x (VMS) b 


Ionization State 


A» (A) c 


f d 


C 


0.1 Mq 


4.1 M 


CII 


1334.5 


0.1278 




C IV 


1548.2 


0.1908 








C IV 


1550.8 


0.09522 





0.5 


44 


O I 


1302.2 


0.04887 


Si 


0.06 


16 


Si II 


1304.4 


0.094 








Si IV 


1393.8 


0.514 








Si IV 


1402.8 


0.2553 


Fe 


0.07 


6.4 


Fe II 


1608.5 


0.058 








Fe II 


2344.2 


0.114 








Fe II 


2382.8 


0.300 



a Metal yield in Scalo IMF case. 
b Metal yield in VMS IMF case. 
c Rest-frame transition wavelength. 
d Oscillator strength. 



TABLE 2 

Line Flux Ratios from minimum to maximum. 



observed time 
IMF 


crossing time 
VMS 


crossing time 
Scalo 


1 day 
VMS 


1 day 
Scalo 


O I/C II 1334.5A 
O I/Fe II 2334.2A 
O I/Fe II 2382. 8A 
O I/Si II 1304.4A 
O I/Fe II 1608. 5A 


0.078-0.086 
0.014-0.030 
0.019-0.031 
0.045-0.059 
0.078-0.090 


0.918-0.965 
0.408-0.444 
0.424-0.458 
0.807-0.881 
0.800-0.852 


0.079-0.090 
0.014-0.042 
0.018-0.045 
0.045-0.089 
0.078-0.105 


0.091-0.097 
0.040-0.043 
0.041-0.045 
0.081-0.088 
0.080-0.085 




VMS 



c 11 
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Fig. 10. — Total GRB spectrum observed at the reverse shock 
crossing time: tgg = 16. 7x (1 + 16.4) s. We show the spectral region 
redshifted into the near IR, accessible to the NIRCam and NIRSpec 
instruments on board the JWST. Metal absorption lines are im- 
printed according to the Pop III SN event, PISN vs. core-collapse 
(Type II). The former originates from a very massive star (VMS) 
progenitor, whereas the latter from a less massive one. In each 
case, the cutoff at short wavelengths is due to Lyman-a scattering 
in the neutral IGM. 

In Figure [TDJ wc show the part of the resulting after- 
glow spectrum that has been redshifted into the near-IR, 
where the JWST is most sensitive. For the two assump- 
tions regarding the Pop III IMF, top-heavy (VMS) and 
more normal (Scalo), the spectra are depicted at the re- 
verse shock crossing time. The metal absorption lines 
are imprinted on the smooth underlying GRB afterglow 
spectrum, rendering their identification straightforward. 
The cutoff at short wavelenghts is due to Lyman-a scat- 



FlG. 11. — Spectral signature of metal absorption lines at the 
reverse shock crossing time: tfg = 16.7 X (1 + 16.4) s. We use the 
metal yields in the VMS case. The observed equivalent widths 
(EWs) are also shown in the figure. If multiple lines are present, 
the given EWs belong to the strongest lines in each panel. 

tering in the IGM which is expected to be still completely 
neutral at z > 10. For detailed predictions, it is ad- 
vantageous to quantify the strength of absorption lines 
through equivalent widths (EWs). The observed EWs 
can be calculated from 

W=(l + z) J [I- e- T ^]dX . (22) 

In Figures [IT] and I12[ we zoom in on the vicinity of 
the strongest lines, again at the reverse shock crossing 
time for both IMF cases. In addition, we explicitly give 
the derived EWs, which are about a few ten A for a 
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Fig. 12. — The same as Fig. 1111 but displaying metal absorption 
lines for a normal (Scalo) Pop III IMF. Individual EWs are again 
given for the strongest lines in each panel. 

top-heavy IMF, and an order of magnitude smaller for a 
normal (Scalo) IMF. The EW values do not change sig- 
nificantly, if evaluated at 1 day after the explosion. Such 
line strengths are well within reach of the NIRSpec in- 
strument on board the JWST, which has a spectral res- 
olution up to R = A/AA = 1000, and a line sensitivity 
of 10 -18 ergs -1 cm~ 2 (Gardner et al. 2006). In general, 
the limiting equivalent width of an unresolved line at ob- 
served wavelength A b s , at 5<7 significance, is (Tumlinson 
et al. 2002): 

w - - WW) ■ < 23 > 

where S/N is the signal-to- noise ratio per resolution el- 
ement. Within our assumptions, therefore, spectroscopy 
of Pop III GRB afterglows will be able to distinguish 
between different enrichment events, PISN vs. core- 
collapse. There is, however, a possible complication. The 
total enrichment level would be rather similar to that 
from one PISN for the case of 10 Type II SNe. Such a 
situation is not implausible, given the number of Pop III 
star formation sites, i.e., minihalos, in the vicinity of the 
Pop III GRB (see Sec. 3). However, we could then still 
discriminate between the cases using absorption line flux 
ratios, which remain quite different. 

In Table 2, we list ranges for select flux ratios, as en- 
countered along different LOSs. We consider the ratios 
O I/X, where X denotes another species, at both the re- 
verse shock crossing time and at 1 day after the explosion. 
In the former case, there is an order of magnitude differ- 
ence between the two IMFs, allowing a robust identifi- 
cation. At later times, however, this clear-cut signature 
disappears again. In Figure Q21 we show the cumulative 
EW distribution for 100 randomly selected sightlines at 
redshift z — 16.4. Again, the resulting distributions are 
quite different for the two IMF cases. It is, however, im- 
portant to note that the statistical significance of this 
analysis is limited by the fact that it is based on only 
one simulation, and only one GRB afterglow. In reality, 
observations will sample multiple explosion sites with re- 
spective LOS distributions that sample many clouds at 
various evolutionary stages along the way. Our analysis 



Fig. 13. — Cumulative distribution of metal line equivalent 
widths. The curves are calculated by tracing 100 randomly selected 
sightlines through the simulation box, centered on the location of 
the GRB. Top panel: VMS case. The curves show the main lines, 
as labelled. Bottom panel: Normal IMF case, corresponding to 
conventional core-collapse (Type II) SN explosions. It is evident 
that the resulting line statistics is quite distinct for the two IMF 
cases. 

should still provide a representative view of the typical 
situation, whereas a more complete study has to await 
the further increase in available computational power. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we develop a diagnostic to study pre- 
galactic metal enrichment in the vicinity of the first 
galaxies. Specifically, we utilize the bright afterglow of a 
Pop III GRB as a featureless background source, and cal- 
culate the strength of metal absorption lines that are im- 
printed by the first heavy elements produced by Pop III 
SNe. To approximately derive the metal absorption line 
statistics, we use an existing highly-resolved simulation 
of the formation of a first galaxy which is characterized 
by the onset of atomic hydrogen cooling in a halo with 
virial temperature > 10 4 K. The reverse shock initially 
dominates the afterglow flux a few hours after explosion, 
followed by the forward shock emission later on. The 
fluxes in the near-IR and radio bands may be detectable 
with the JWST and the VLA out too 30. We pre- 
dict that the afterglow emission peaks from near-IR to 
millimeter bands with peak fluxes from mJy to Jy at dif- 
ferent observed times. Metal absorption lines in the GRB 
afterglow spectrum, giving rise to EWs of a few tens of 
A, may allow us to distinguish whether the first heavy 
elements were produced in a Pop III star that died as 
a PISN, or a core-collapse SN. To this extent, the spec- 
trum needs to be obtained sufficiently early, within the 
first few hours after the trigger. 

The absorption signature of the first SN events might 
allow us to constrain the underlying mass scale of Pop III 
stars. This is the key input parameter to predict their 
evolution, nucleosynthetic yields, and modes of death, 
which in turn determine how the first stars impact sub- 
sequent cosmic history. The strength of the associated 
stellar feedback governs the assembly process of the first 
galaxies, in the sense that the stronger feedback from 
more massive stars shifts their formation to later stages 
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in the hierarchical build-up of structure (e.g., Ricotti et 
al. 2002; Greif et al. 2010; Frebel & Bromm 2012; Ritter 
et al. 2012; Wise et al. 2012). 

It will, however, be very challenging to directly probe 
the initial epoch of cosmic star formation. The reason 
is that even the JWST will not be able to detect indi- 
vidual first stars, but instead only more massive stellar 
systems or clusters that form later on in more massive 
host systems (e.g., Pawlik et al. 2011). Pop III GRBs 
may thus afford us one of the few direct windows into 
the crucial epoch of first light, another one being the 
extremely energetic pair-instability SN or hypernova ex- 
plosions that are predicted to end the lives of the most 
massive Pop III stars (e.g., Pan et al. 2012; Hummel et 
al. 2012). The promise provided by the combination of a 
wide field GRB trigger mission, such as JANUS or Lob- 
ster, with the next-generation of highly-sensitive near- 



IR telescopes, such as the JWST or the planned ground- 
based extremely large facilities (the E-ELT0, GMT0, and 
TMT^I), is huge. GRBs are likely to play a key role in 
finally opening up the high-redshift frontier, all the way 
back to the very beginning of star and black hole forma- 
tion. 
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